Phase diagram for condensation of microcavity polaritons: from theory to practice 
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The first realization of a polariton condensate was recently achieved in a CdTe microcavity 
[Kasprzak et al, Nature 443, 409 (2006)]. We compare the experimental phase boundaries, for 
various detunings and cryostat temperatures, with those found theoretically from a model which 
accounts for features of microcavity polaritons such as reduced dimensionality, internal composite 
structure, disorder in the quantum wells, polariton-polariton interactions, and finite lifetime. 
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I. INTRODUCTION 



Microcavity polaritonsii^ are the quasi-particles which 
result from strong light-matter coupling in semiconductor 
microcavities. The light polariton mass means a weakly 
interacting dilute Bose gas picture would imply a high 
transition temperature^i^i^ and, even with more sophis- 
ticated methods,— iii^ the expected transition tempera- 
ture is 10^ times higher than for atomic gases. While 
Bose-Einstein condensates (BEC) of atomic gases are 
now routinely studied, condensation of microcavity po- 
laritons has been long sough1^ i^°'^^i^^'^'^i^^i*^ and only re- 
cently has a thermal equilibrium condensate accompa- 
nied by spontaneous coherence been realized in a CdTe 
microcavity In the same structure the polarisation has 
been investigated;^ second order temporal coherence has 
been studiedfi^ and pinned vortices have been observed 
in the condensed phase.— Spontaneous coherence of mi- 
crocavity polaritons has been also recently investigated 
in GaAsr^i^i while there have been early reports on room 
temperature polariton lasing in a GaN cavity 

Microcavity polaritons are interacting particles in a 
disordered potential with internal structure and finite 
lifetime, which form a finite two-dimensional gas, so their 
condensation differs from that of ideal non-interacting 
three-dimensional (3D) bosons for which BEC was orig- 
inally discussed. It is therefore crucial to establish these 
differences and the nature of polariton condensation. To 
this end, the aim of this paper is to compare theoretical 
and experimental results for the phase diagram, and to 
discuss what they reveal about the nature of the polariton 
condensate. We discuss both complete thermal equilib- 
rium phase diagrams, and also consider effects that may 
result from pumping and decay (even when the polari- 
ton distribution is well thermalised). As a further way 
of illustrating the differences between the ideal 3D Bose 
gas and microcavity polaritons, the condensate fraction 
as a function of pumping strength is presented; this is 
compared qualitatively to a theoretical model including 
effects of pumping and decay. 

We experimentally evaluate the critical density for 



condensation at low cryostat temperatures in the range 
Tcryo = (5, 15)K and for exciton-photon detunings in the 
range S = (5.3, 12.5) meV. These densities arc compared 
with the predictions of a realistic microscopic model that 
accounts for the composite structure of microcavity po- 
laritons, interactions enhanced by quantum well disorder 
and finite polariton lifetime. This comparison reveals 
two crucial aspects of the current realisation of polariton 
condensation in CdTe: that it lies in the crossover be- 
tween BEC of weakly interacting bosons and the mean- 
field phase transition driven by interactions; and that the 
excitonic disorder is important in the evaluation of the 
polariton-polariton interactions. 

The organisation of the remainder of this paper is as 
follows: Section [ll] describes the experimental estimate 
of critical density, found by integrating the occupation. 
Section|TTT]then introduces the theoretical models we con- 
sider, and in Sec. IIII CI we discuss the various theoreti- 
cal estimates of phase boundaries they produce. The 
comparison between the various estimates of experimen- 
tal critical density and the theoretical phase boundary is 
shown in Sec. IIVI The effects of pumping and decay are 
discussed in Sec.|Vl and their influence on the condensate 
fraction - in Sec. IV Al 



II. EXPERIMENTAL DETERMINATION OF 
CRITICAL DENSITY 

A CdTe/CdMgTe microcavity with 16 quantum wells 
and Rabi splitting f2fl = 26meV, is non-resonantly 
pumped at the first high energy lobe of the Bragg struc- 
ture, at 1.768eV (about ^ lOOmeV above the lower 
polariton ground state energy) by a continuous-wave 
Ti:sapphire laser combined with an acousto-optic modu- 
lator (1/xs pulse duration with a duty cycle of 1%) to re- 
duce the sample heating. The excitation spo t is uniform 
and of ~ 35//m diameter. Following Ref[la, we detune 
the photon energy Uk above the exciton iJx [S = loq — E-^) 
and, increasing the pump laser power, we measure the 
energy- and angle- (or momentum-) resolved emission in- 
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tensities via far-field spectroscopy (inset of Fig. [T]). In- 
plane momentum k and emission angle 9]^ are related by 
Hc\h\ ~ i?^^ sin 0k, where 



LP,UP 



are the lower (LP) and upper polariton (UP) energies. 
From the energy-integrated photoluminescence (PL) at 
a given momentum T^k, the occupation Ok of that state 
is given by^ 



cos"* 6'kC^ ' 



Ck 



1 



t^k - -Ex 



By finding the energy at which, for a given momen- 
tum, the PL is maximum, momentum dependence can 
be converted to energy dependence; an example is shown 
in Fig. [T] Above a threshold pump power, the ground 
state occupation grows exponentially, while the effec- 
tive LP temperature, Tpff, (extracted from the tail of 
the occupation) docs not change much. This macro- 
scopic occupation of the ground state has been shown 
to be associated with macroscopic coherence across the 
spot size,— demonstrating polariton condensation. At 
threshold T^s = 16K; this is above Tcryo because the po- 
lariton lifetime (of l-2ps) is shorter than typical relax- 
ation times (although longer than their thermalization 
time)J^ii^ The ground state occupation at threshold is, 
within a factor of two, measured to be one; in Fig. [l]this 
is rescaled to one. 
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FIG. 1: (Color online) LP occupation versus energy at various 
excitation powers extracted from the energy-resolved far-field 
emission (inset). 

The LP density at the threshold pump power can then 
be estimated by integrating the occupation to give the 



total density. The estimate found this way is represented 
on Fig. [3] by blue stars. Since the system is thermalized, 
by rescaling the occupation at k = to one, this density 
estimate is determined by the effective temperature only; 
any variation at fixed temperature is due to changes to 
the LP density of states. 



III. THEORETICAL PHASE DIAGRAMS 

This section discusses the theoretical modelling of the 
phase diagram, and further methods of estimating the 
experimentally observed critical density. The discussion 
here concerns the thermal equilibrium case; possible ef- 
fects of pumping and decay are discussed later in Sec. |Vl 
The discussion will compare the use of two models; the 
first — model ([T]) — is a model of disorder localised sat- 
urable excitons coupled to propagating photons.— At 
low enough densities and temperatures, the critical den- 
sity given by model ^ will be dominated by long wave- 
length bosonic fluctuations^''^ and so the results of this 
model can be reproduced by modelling lower polaritons 
as a weakly interacting dilute Bose gaS j^'^'^i^^ which we 
shall call model ([2]). Excitonic disorder can significantly 
modify the effective polariton interaction compared to 
the clean case. To take account of excitonic disorder in 
model ©, we choose the interaction strength for this 
model so as to match the low density limit of model ((T|) , 
in which excitonic disorder is explicitly considered. Such 
matching is achieved by calculating the form of blueshift 
vs density in both models, and choosing the effective in- 
teraction strength for the bosonic model so they agree. 
This calculation of blueshift vs density also allows a sec- 
ond estimate of the experimentally found critical density, 
by measuring the blueshift observed at that density. 



Localised exciton model 



The first model, following Refsl7ll24l. starts from the 
single-particle exciton states (numerically) evaluated in 
a disordered quantum well, as explained further below. 
Although being 2D, all such states are localised, the na- 
ture of the exciton states and their associated oscillator 
strength change substantially with the exciton energy i21 
The exciton states that couple most strongly to the cav- 
ity photons are those just below the band edge. 

Although the exciton states are localised, polaritons 
consist of a superposition of many excitons^^ and so need 
not be localiscdj2£ Polaritons in CdTc are also localised 
by photonic disorder, associated with the spatial inho- 
mogeneity observed in the PL above threshold. However, 
above threshold, coherence can be observed over a large 
length scale approaching the excitation spot size^^i^ and 
so polaritons could not be localised on a shorter length 
scale. This gives an energy scale up to 0.1mcV~ IK, 
much smaller than the observed blue-shift (due to inter- 
action). This difference in energy scales suggests that 
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localisation is a perturbation to a theory of interacting 
bosons, rather than vice versa; we therefore neglect pho- 
tonic disorder when finding the phase boundary. 



structure of the exciton smooths disorder on lengthscales 
shorter than the exciton Bohr radius^l Ox, so we can con- 
sider a Gaussian noise correlated on a length scale Ic'- 



1. Many body Hamiltonian 

In the many body Hamiltonian, we approximate the 
interactions between excitons by exclusion: Each single- 
particle exciton state can only be occupied once. This 
approximation over-estimates the on-site interaction and 
neglects the inter-site interaction — in reality there is 
finite Coulomb and exchange energy to multiply occupy 
the same single exciton states. However, as discussed in 
Refii, this approximation is valid at low enough densi- 
ties. This procedure yields model ([1]), 

a k 

+ 4f 5^ (•9"'kV'k^+ + h.c.) . (1) 



(F(R))=0, (y(R)l/(R')) = 



A 



iq-(R-R') 



The spin operators S^,S^ describe states |ta) and \la) 
corresponding respectively to the presence or absence of 
an exciton of energy . The operators ip^. describe con- 
fined photon modes of energy Wk = -t- ?i^k^/2mph and 
Wph/wo = 2.58 X 10~^ (mo is the bare electron mass). 
The quantisation area A that appears in Eq. ^ plays 
no role in any final answer, and is taken to infinity in 
all our calculations. The Hamiltonian written in Eq. ([1]) 
does not yet take account of spin of the exciton states or 
polarisation states of the light. To a first approximation, 
the interaction between different spins is weak; in such a 
limit each polarisation behaves separately, and one has 
independent copies of this model for each spin, thus in- 
troducing a degeneracy factor = 2 in the total density. 
Beyond this limit, opposite spin polarisations are weakly 
attractive, producing bound excitonic states; the effect 
of such an interaction is discussed in Appendix [Xl 



2. Energies and coupling strength of excitons 



As in Refsl7ll24l. the exciton energies Sa and coupling 
strengths ga^k (which lead to the inhomogencous exci- 
ton lincwidth) are taken from numerical diagonalisation 
of exciton wavefunctions in a disorder potential^i These 
numerical calculations are described in detail in 
so we only brieffy summarise the method here. We nu- 
merically solve the Schrodinger equation in a disorder 
potential: 



.3. 

2 m,, 



y(R) + 



$„(R) =e„$a(R) 



We use the CdTe exciton mass^ m^/mo = 0.6, and 
a Gaussian correlated exciton disorder. For disorder 
weaker than the exciton binding energy, the internal 



We take £c = 167A > a-^ and so the only free parameter 
in this calculation is the strength a of the disorder; this 
was chosen to give an exciton linewidth that matches the 
observed CdTe exciton lincwidth {^-^ ImeV), leading to 
cr = 0.79meV. 

The exciton-photon coupling strength involves the 
overlap between the disorder localised exciton state, and 
the photon wavefunction. Since the relevant photon mo- 
menta are small compared to the typical scales of the 
exciton wavefunction (i.e. small compared to 1/ax), we 
take this coupling strength at zero photon momentum, 
as ga,k — ga.o (*&Q|k = 0), where $q, is the disorder 
localised exciton state. To fix the overall scale of the 
couplings we use Rabi splitting fii? = 2^yj2a \9a,o\'^/A 
with the observed Rabi splitting. Thus, the complete set 
parameters defining model ([1]) are the exciton mass, the 
Rabi splitting, the exciton linewidth, the exciton Bohr 
radius and the photon mass. These parameters all have 
reasonably well established values, and arc not used as 
fitting parameters. 



B. Weakly interacting dilute Bose gas 

At low densities and temperatures, the thermodynam- 
ics of model (HJ is dominated by long wavelength bosonic 
fiuctuations,— and so it gives the same predictions as 
an effective LP model^'^^'^^ [model 



k,k',q 



(2) 

As before, in the absence of biexcitonic effects, the spin 
degree of freedom only contributes a factor gs ~ 2 
to the density. Such a model for the condensation of 
weakly interacting bosons is often studied in the con- 
text of excitons, and of atomic gasesi^ For polaritons, 
model ^ has been considered^ii with the interaction 
strength V^^^^., ^ taken from a combination of Coulomb 
interaction and the nonlinearity of exciton-photon cou- 
plingi^ However, excitonic disorder can significantly en- 
hance polariton-polariton intcractionsii^ and make them 
larger than those in clean QWs. In choosing the value 
of V^^ij., q we wish to include these effects of disorder, 
which are described microscopically in model ([T]). As 
mentioned above, the form of model ^ describes the 
relevant excitations of model ([T]) at low densities and 
temperatures, and so the physical predictions of these 
two models should match in this limit. Hence, the effect 
of disorder on the interaction strength can be included 
by choosing the value of V^'^., ^ so as to make the results 
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of model ([T]) and model ([2|) mateh at low densities and 
temperatures. 



1. Fixing effective interaction strength 

Following the relation between the two models as de- 
scribed above, the value of V^^ is chosen so that the 
blueshift vs density in the normal state calculated in both 
models will match. In the normal state of model ^ , the 
LP blue-shift varies linearly with the density 5E^^ ~ 
V^^n. In model ([1]) the blue-shift depends on both tem- 
perature and disorder, but, in the low density normal 
state, it also varies approximately linearly (lower panel 
of Fig. [2]). We choose so these gradients are the same. 
In both models, once condensed, the LP branch locks to 
the chemical potential, and so the blue-shift traces the 
dependence of chemical potential on density, increasing 
more rapidly then linearly at first, before eventually sat- 
urating (illustrated by the red x symbols in the lower 
panel of Fig. [2]). 
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FIG. 2: (Color online) Comparison of (panel a) experimental 
LP blue-shift (red data with error bars) vs. incident pump 
power; (bottom panel) theoretical LP blue-shift (red x sym- 
bols) vs. mean-field density. The experimental LP linewidth 
(black dashed) and threshold (found from the occupation 
data) are also shown. Measured blue-shifts at threshold can 
be translated to densities as explained in the text. 

The value found by this procedure is V^^ = 
1.52/xeV/xm-^, implying a blue-shift of 0.304meV at a den- 
sity of 10^°cm~^). This effective interaction strength V^^ 
is roughly 10 times that found from Coulomb and satu- 
ration effects in a clean systemi^ The relative size of 
inter-site Coulomb and saturation effects can also be seen 



in the UP energy shift with pump power: The Coulomb 
term shifts both LP and UP in the same direction, satura- 
tion effects lead to opposite shifts. Experimental data in 
the on-line supplementary information of Ref[l^ show an 
UP red-shift of magnitude comparable to the LP blue- 
shift, strongly suggesting that saturation effects domi- 
nate Coulomb interactions. 



2. Estimating experimental density from blueshift 

The calculation of blue shift as a function of density 
allows one to extract a second estimate of the LP den- 
sity at threshold, by converting the experimentally ob- 
served blueshift (Fig. [2l upper panel) to the correspond- 
ing density. However, as seen in Fig. [21 the mean-field 
estimate of blue-shift differs according to whether one 
considers the condensed or uncondensed case. Since fluc- 
tuations beyond the mean-field approximation delay the 
transition, the actual blue-shift at threshold should lie 
between the condensed and uncondensed mean-field es- 
timates. Thus, for a given observed blue shift, one can 
extract two estimates of the density that should bound 
the actual density. These two translations to density are 
shown in Fig. [31 as red upward triangles for the con- 
densed calculation, and green downward triangles for the 
uncondensed case. 



C. Theoretical phase boundaries 

The experimental estimate of the critical polariton 
density can be now compared to the theoretical phase 
diagram. From model ([T]), one may derive a mean- field 
phase diagram, for condensation of polaritons in the low- 
est energy state. The mean-field phase boundar y 
comes from finding the mean-field estimate of density 
(gs = 2 takes into account the degeneracy introduced by 
the excitonic spin): 



-y 

A ^ 



A 



1 — tanh 



(3) 



at the chemical potential \i which first satisfies the saddle- 
point equation (which plays a role analogous to the self 
consistency condition or Gross-Pitacvskii equation for a 
weakly interacting dilute Bose gas): 



-y 



A 



.9a,0| 



;tanh[/3(e - ^)/2] 
e — [I 



(4) 



This mean-field boundary is shown as a black solid line 
in Fig. [3 

Including fiuctuation corrections to this mean-field 
transition, one recovers a smooth crossover from the 
mean-field limit at higher densities to a fiuctuation domi- 
nated regime at low densities^i^; in this fiuctuation dom- 
inated regime, the results of model (|T|) are equivalent to 
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those of model ([2]). This crossover is at a temperature 
that is controlled by Rabi splitting, and is of the order 
of one tenth of the Rabi splitting. In 2D, for interact- 
ing bosons with a quadratic dispersion, the transition 
temperature to the superfluid state, the Berezhinskii- 
Kosterlitz-Thouless (BKT)^^ temperature, varies lin- 
early with the density n and is close to the quantum 
degeneracy temperature, Tdeg oc n/m. The results of 
model dl]) are not quite so simple as the dispersion is 
non-quadratici^i^ However, at ultra-low densities, the 
phase boundary of model ^ (blue dotted-dashed line) 
recovers a linear dependence of critical temperature on 
density (black dashed Hne). 

The BKT transition discussed above describes an in- 
finite two-dimensional system. The experimental system 
however contains photonic disorder, which might act as 
a trapping potential for polaritons. In a non-interacting 
harmonically trapped two-dimensional Bose gas, a tran- 
sition exactly like that described by Bose and Einstein 
can exist due to the modification of density of states by 
the harmonic trap^i^ However, the polariton system is 
interacting, and it is not a-priori valid to ignore inter- 
actions; the inclusion of interactions in an harmonically 
trapped system^! can replace the BEC transition of the 
non-interacting gas with a transition better described by 
BKT physics. The exact details of the critical behaviour 
in a trapped two-dimensional interacting system are del- 
icate. As already discussed in Sec. IIII A| for the parame- 
ters relevant in the current system, the effects of polari- 
ton trapping induced by the photonic disorder are weak 
compared to the polariton-polariton interaction strength. 
However our interest here is in the prediction of the phase 
boundary, we note that for the parameters relevant here 
simple estimates of the critical temperature for the BEC 
of non-interacting trapped system do not differ signifi- 
cantly from those for the BKT transition. There is a 
quite separate question concerning coherence (and con- 
densate fraction) in a two dimensional trapped system, 
which we discuss below in Sec. IV Al 




n [cm ] n [cm ] 



FIG. 3: (Color online) Theoretical phase diagram of temper- 
ature versus density (gray region is condensed) for different 
values of the detuning S and Tcryo, with superimposed the 
estimates of density from the PL (blue stars) and from the 
blue-shift of LP (red upper and green lower triangles) . 

rather than with the condensation transition itself. As 
an upper bound on this source of uncertainty wc note 
that the density when the lincwidth reaches a minimum 
is between 1.6 and 2 times that at threshold. For the 
estimates from the blueshift, the uncertainties here are 
harder to quantify. The exclusion interaction for excitons 
may overestimate blueshift, implying our estimated den- 
sity is too low. Alternatively, our calculation assumes a 
thermalised polariton distribution, whereas the influence 
of the exciton reservoir and hot carriers remains unclear 
within the present modelling, these might increase the 
blueshift, suggesting that our estimated density is too 
high. Finally, the theoretical phase boundary as yet ne- 
glects dephasing; it is this effect to which we next turn. 



V. EFFECTS OF PUMP AND DECAY 



IV. COMPARISON OF THEORETICAL AND 
EXPERIMENTAL BOUNDARIES 

In Fig. [3] the experimental density is near the crossover 
between the low- and high-density parts of the theoret- 
ical phase boundary. The density estimated from the 
blue shift is slightly greater than that from integrating 
the occupation; however without the effects of excitonic 
disorder, the density from the blueshift would be ten 
times higher. Let us consider the possible systematic 
errors involved in estimating the density. In integrating 
occupation, the sources of uncertainty are the normali- 
sation (discussed earlier), and the accuracy with which 
the critical pumping power has been determined. We 
may underestimate the experimental critical density if 
the nonlinear threshold we use is in fact associated with 
the diverging susceptibility on approaching the transition 



As is apparent from Fig. [H and the discussion in 
Reffl6l. the polariton distribution is thermalised. How- 
ever a thermalised population distribution docs not nec- 
essarily imply that one can neglect effect a"'^^i'^^'^^i'^°'^"'^i"^^ 
of pumping and decay. To understand this distinction, 
let us consider two dimensionless ratios that characterise 
non-equilibrium effects. One might estimate the effect 
of pumping and decay by considering the ratio of po- 
lariton lifetime to thermalization time (e.g., ~ 2 in the 
current experiment, ^ 2000 in early atomic gas^ BECs). 
When this ratio is larger than one, the distribution will 
be thermalised. However, even with a thermal polariton 
distribution, pumping and decay cause a particle flux, 
introducing dephasing, which increases the critical den- 
sity at a given temperature.— The importance of this 
dephasing can be estimated by comparing the homoge- 
neous photon linewidth, LW~ ImeV, to the tempera- 
ture ksT w 2meV. As these parameters are comparable. 
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lincwidths will have significant homogeneous contribu- 
tions, and dephasing can affect the ability of the system 
to maintain a coherent state. 

As an estimate of how such dephasing would mod- 
ify critical density at a given temperature, Fig. [4] com- 
pares the phase boundaries found from calculations using 
model ([1]), but with photon decay and injection of exci- 
tons^i^S taken into account. These calculations intro- 
duce a flux of particles by allowing photons to decay and 
continually injecting new excitons to maintain a steady 
state. The departure from the closed system picture 
is controlled by two parameters, the photon decay rate 
2k = ImeV, and a pumping rate 7 that describes cou- 
pling between the polariton system and the reservoir of 
high energy excitons. The value of 7 is not clearly known, 
but is bounded. For 7 3> k one has a polariton linewidth 
larger than the bare photon linewidth; for 7 <C k there is 
instability of the homogeneous condensate^; a range of 7 
bounded by these consideration is shown for illustration 
in Fig. |4] — within these limits, the exact value of 7 does 
not significantly shift the curves. Although these calcula- 
tions use model ([T} , for tractability of the nonequilibrium 
calculations, a Gaussian distribution of excitonic energies 
£a and a constant coupling ga,-k = are used. The 

main effect of this replacement is that there is no simple 
relation between the units of density for this calculation 
and those in the previous figures, hence the arbitrary 
units in Fig. |4l However, one may note that the form of 
the equilibrium boundary (black) in Fig. 0] is very simi- 
lar to that in Fig. [31 and one may thus estimate the size 
of shift associated with pumping and decay. Dephasing 
shifts the phase boundary to higher densities; this would 
reduce the distance between the predicted phase bound- 
ary and the critical density estimated from the blueshift. 



A. Condensate fraction 

As a further illustration of how dephasing depletes the 
condensate, we compare experimental condensate frac- 
tion (Fig. [S]) with the theoretical quasi-condcnsate frac- 
tion (evaluated at the mean-field level) with and without 
pumping and decay (Fig. [4]). Although trapping negli- 
gably affects the phase boundary (the properties of the 
phase boundary being dominated by interactions), trap- 
ping and finite size are important in allowing a non-zero 
condensate fraction4^i^ A quasi-condensate fraction can 
also be defined for an infinite two-dimensional system; its 
definition follows from the behaviour of the off-diagonal 
one-particle density matrix p{r). In three dimensions 
p{r) tends to a constant at long distances, which defines 
Off-Diagonal Long-Range order; in two dimensions it in- 
stead decays as a power law. However, the power law 
decay is due to phase fluctuations, and by separating the 
effects of phase fluctuations from the quasi-condensate 
depletion due to density fluctuations^ one may extract a 
quasi-condensate density. In a trapped two-dimensional 
system there can be a regime near the phase transition 
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FIG. 4: (Color online) Critical temperature versus density 
(in arbitrary units) (panel a) and quasi-condensate fraction 
at T = 20K (panel b) including the effects of pumping (rate 
7) and decay, with decay rate for a homogeneous cavity line- 
width of ImeV, for a model Gaussian distribution of the ex- 
citonic energies. In panel (a), for 7 = 0.97meV and high 
temperatures the nonequilibrium model has an instability of 
the homogeneous condensate^; the boundary of this instabil- 
ity is marked by the dashed line, and the stable condensed 
region is marked by hatchings. 

in which there is quasi-condensation; however as tem- 
perature is reduced or density increased, coherence will 
rapidly extend over the entire systemi^i^ Thus, apart 
from a small region near the transition, coherence across 
the entire system is expected even when the phase tran- 
sition itself is like the BKT transition. 

The experimental condensate fraction is determined 
from the occupancies illustrated in the inset of Fig. [5] 
as follows: The occupancy shows a clear peak at low mo- 
mentum; the experimental estimate of condensate frac- 
tion comes from the fraction of polaritons found in this 
low momentum peak, i.e. the integrated density within 
the FWHM of the occupancy peak, see Refl47l for further 
details. The theoretical mean-fleld estimate is the ratio 
of the coherent density, + J2a K-^q )P *° total 

(mean-fleld) density of polaritons 1^ There are of course 
complications in comparing the experimental condensate 
fraction of a trapped 2D Bose gas^^ to the theoretical 
mean-fleld condensate fraction; however the comparison 
between the various mean-field estimates clearly shows 
that pumping and decay arc responsible for a reduction 
of the quasi-condensate fraction. 

VI. CONCLUSIONS 

Concluding, we have compared experimental and the- 
oretical phase boundaries for condensation of CdTe mi- 
crocavity polaritons. The experimental data lie near the 
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FIG. 5: (Color online) Experimental condensate fraction (the 
fraction of polaritons in the low momentum peak), at de- 
tuning 5 = 10.5meV, Tc^yo = 5.4K, T^s = 16.5K. Inset: 
Occupancy vs angle, showing the peak at small angles (low 
momentum), the integrated weight inside which is taken as 
condensate density. 



crossover between a regime where the transition can ad- 
equately be described by condensation of structureless 
bosons and a mean-field regime where instead the long- 
range nature of the polariton-polariton interaction deter- 
mines the boundary. We were limited to temperatures 
near this crossover region: Because of the short polari- 
ton lifetime, the polariton temperature is decoupled from 
that of the lattice. There are small differences between 
the various ways of estimating experimental density; were 
one to neglect effects of disorder these differences would 
have been an order of magnitude larger. However, these 
estimates of experimental density lie close to the theoret- 
ical phase boundary, and considering the effects of pump 
and decay improves this agreement. 
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APPENDIX A: INTERACTION BETWEEN 
SPINS AND BIEXCITONS 

This appendix discusses the effect on the phase dia- 
gram of including an attractive interaction between op- 
posite polarisations of two-level systems. The model used 
is a straightforward generalisation of Eq. ([1]), explicitly 
keeping track of the four possible configurations of each 
disorder-localised site; unoccupied, singly occupied by an 
exciton or occupied by a biexciton, thus each site has an 
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FIG. 6: (Color online) Mean-field phase boundary calculated 
with varying strengths of biexciton binding energy, with a 
detuning 5 — 5.3meV (i.e. corresponding to top-left panel of 
Fig.0. 



effective Hamiltonian: 



/ 

a;. 



A. 



A 



a,R 




V 



a,R 









A. 
A„ 



A* 



A 



a,R 



where A 



{L,R), 



'■a,R 

(Al) 

Ek5a,ki/'(L,fl)k/V^ describes the cou- 
pling of a given localised site to the left- or right- 
circularly polarised photon field, and -Bxx is the biexciton 
binding energy. This expression replaces the combination 
[Sa - + I]k {9ot.yi'\^S'^ + h.c.) / that appears in 
Eq. H]) for each site. The mean-field theory, analogous to 
Sec, nil Cl is found by calculating the mean-field estimate 
of density which may be written as: 



2(za zlX) 



(A2) 



where we have defined: 



A 



and the gap equation (found from the linear polarisability 
of the four-level system) becomes: 



Wo - M 



\- 2Za + Zl\ 



1 



- [I - Exx 

(A3) 

It is straightforward to see that if -Exx = 0, these two 
forms reproduce those in Eq. ([3]) and Eq. @ respectively. 
A more complete treatment of binding ought to take ac- 
count of how the binding strength varies according to the 
localised nature of the single exciton states, however to 
place a bound on how large a shift of critical tempera- 
ture the biexciton binding might induce it is sufficient to 
consider this model with a constant biexciton binding. 

Figure [S] shows the effect with binding energies up to 
8meV, corresponding to the largest value seen for exci- 
tons bound at vicinal surfaces in CdTe/CdMgTe wells4^ 
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Such a value is almost certainly an overestimate for the 
samples used in current experiment, as the results of 
Refl48l are for narrower quantum wells than the current 
experiment (3.2nm vs 5nm), and correspond to biexci- 
tons with the greatest in-plane confinement. Even with 
this overestimate of biexciton binding, the effect of biex- 
citon binding on the critical temperature can be seen 
to be small. Moreover, were biexcitons to play a sig- 



nificant role, one would expect to see a crossover be- 
tween the photoluminescence associated with exciton- 
polaritons and that from biexciton-polaritons, as was in- 
vestigated in GaAs4^ The fact that no such crossover 
is seen, and that the effects of biexciton binding on the 
mean-field phase boundary are small justify our neglect 
of biexciton physics in the rest of the paper. 
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